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We use recently developed strong-coupling expansion methods to study the two-particle spec- 
£f>^ | tra for the frustrated alternating Heisenberg model, consisting of an alternating nearest neighbor 

■ antiferromagnetic exchange and a uniform second neighbor antiferromagnetic exchange. Starting 
from the limit of weakly coupled dimers, we develop high order series expansions for the effective 
Hamiltonian in the two-particle subspace. In the limit of a strong applied dimerization, we calculate 
accurately various properties of singlet and triplet bound states and quintet antibound states. We 
also develop series expansions for bound state energies in various sectors, which can be extrapolated 
using standard methods to cases where the external bond-alternation goes to zero. We study the 
properties of singlet and triplet bound states in the latter limit and suggest a crucial role for the 
bound states in the unbinding of triplets and deconfinement of spin-half excitations. 

PACS Indices: 75.40Gb, 75.10Jm, 75.50Ee 
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Q ; I. INTRODUCTION 

In recent years the subject of fractional excitations, or excitations with fractional quantum numbers compared to the 
^ ■ non-interacting limit, has attracted considerable attention. Many experiments in high temperature superconducting 
materials have been interpreted in these terms, and many theories of high temperature superconductors are built 

■ around such fractional excitations. However, fractional excitations in a closed system can only arise in groups that 
have the full quantum numbers of the non-interacting limit. Thus, their unambiguous identification in numerical 

y—i • calculations remains difficult. 

In this work we are interested in studying models, where as parameters in the Hamiltonian are varied, one goes from 
a phase where the excitations have normal quantum numbers to one where they have fractional quantum numbers. 
We would like to develop series expansion methods by which such transitions can be studied and the onset of fractional 
excitations can be demonstrated. Although we work with ID models, the basic methods we develop can be applied 
in higher dimensions as well. 

*^ ■ The best known example of a fractional excitation is a spinon in the spin- half antiferromagnetic Heisenberg chain [W . 
It is well-known that the low-lying excitations of this Bethe-ansatz solvable model consist of a two-spinon continuum. 
Another simple example of a spin-half excitation is a soliton in the Majumdar- Ghosh model It is a domain- wall 
which interpolates between the two dimcrized ground states of the model. The low-lying excitations, for a closed 
system, in this case also consist of a soliton-antisoliton continuum ||. 

Adding a bond-alternation to the exchange constants of the model leads to confinement of the spin-half excitations. 
The elementary excitations now become triplets and the spinons or solitons are bound into pairs. In the unfrustrated 
case, the role of bond-alternation has been studied by mapping onto a massive Thirring model Q). The frustrated 
case, where the ground-states are spontaneously dimerized, has been of considerable recent theoretical interest due 
to its relevance to spin-Peierls systems such as CuGe03. An explicit bond- alternating term in the Hamiltonian can 
be motivated as a mean-field representation of the inter-chain elastic couplings Q . Uhrig et al. ||] and Affleck and 
collaborators have studied the confinement transition for the soliton-antisoliton pairs when such a term is added 
to the Hamiltonian. 

Here, we approach these transitions from the opposite direction. We study these systems in a strong-coupling 
perturbation theory, which begins with the limit of decoupled spin-dimers and treats the inter-dimer couplings as 
a perturbation. In the limit of weakly coupled dimers, the elementary excitations are triplets, which are weakly 
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dispersive. In this limit our strong coupling theory is highly accurate and we can find all details of various two- 
particle bound (and antibound) states. The overall 2-particle spectrum is much richer than that obtained in previous 
studies. Several singlet and triplet bound states and quintet antibound states are found. The number of bound states 
depends on the coupling constants as well as the wavevector. We study the binding energy and the coherence length 
associated with the bound states. We also study the singularity at the critical wavevectors where the binding energy 
goes to zero and the state merges into the continuum. 

Using high order series expansions and extrapolation methods we also study the uniform limit, where the bond- 
alternation term in the Hamiltonian goes to zero. Thus we approach the limit where the triplet excitations break 
up and spin-half excitations become deconfined. Series expansion results show that as the bond-alternation term 
goes to zero in the Hamiltonian, the spectral weights associated with triplet quasiparticles go to zero and the lowest 
lying singlet and triplet excitations become degenerate. These phenomena provide a remarkably clear and simple 
confirmation of the existence of free spin- half excitations in this limit. 

In the frustrated system, the reorganization of the many-body spectra as the system undergoes the deconfinement 
transition presents an interesting puzzle. Since there is a gap, A, to triplet excitations, the two-triplet continuum 
begins at 2A. Thus, in the confined phase, this continuum is separated from the elementary triplet by a second 
gap. However, when the spin-half excitations are liberated, the resulting low energy spectrum consists of a soliton- 
antisoliton continuum, which does not have such a second gap. The consistency of the two pictures requires that in the 
confined phase, between the elementary triplet and the two-triplet continuum there must be a large number of states, 
which upon deconfinement turn into the continuum. On general grounds, these states must include 2, 3, 4, . . ., triplet 
bound states, i.e. states involving an arbitary number of triplets, which must correspond to a soliton-anitisoliton pair 
with arbitrary separation. We discuss insights from studies of two-particle bound states on this issue. 

Another interesting puzzle lies in the spectrum of the Majumdar-Ghosh model (5 = 0, a = 1/2) near k = tt/2. For 
both S = and S = 1, previous studies |||||| have emphasized a bound state below the soliton-antisoliton (s — s) 
continuum. To our knowledge, it has not been noticed that the two triplet continuum (the s — s — s — s continuum) 
falls below the soliton-antisoliton continuum around this wavevector. In our numerical study, we find that except for 
a tiny region very near k = 7r/2, the latter continuum also falls below the bound states. This raises questions about 
the stability of the bound states away from k = tt/2. 

Another puzzle in our studies is how the energy levels might cross each other. On general grounds, one might expect 
that the levels for n-particles with varying n cross each other as the parameters are varied. For small values of the 
perturbation parameter, the energies are arranged in order of increasing n. However, as one approaches A = 1, low 
energy states from each n-sector may appear even below the two-particle continuum. Since n is not a good quantum 
number, it is not clear how this will reflect itself in our perturbation theory. This deserves further attention. 

The organization of the paper is as follows. In section II, we describe the Hamiltonian studied and the various 
parametrizations used. In section III, we study the regime of forced or externally imposed dimerization. This is a 
regime where our series expansions are convergent and we present spectra, binding energies, coherence lengths etc. 
in great detail. In section IV, we consider the regime of spontaneous dimerization, which requires the use of series 
extrapolation methods. In section V, we present discussions and conclusions. 



II. HAMILTONIAN 

We wish to study the alternating Heisenberg chain with frustration ^,^|Jl^-|lJ,^] , 

H = 12 K 1 + (- 1 )^) S » ■ S *+i + aS * • S *+2] . (!) 

i 

where the Si are spin-^ operators at site i, a parameterizes a next-nearest neighbor coupling and S is the alternating 
dimerization. We rewrite the Hamiltonian as 

H = (1 + S) [S M ■ S 2i+ i + A(S 2i • S 2i _! + yS t ■ S l+2 )] , (2) 

i 

The parameter space (5, a) is equivalent to the parameter space (A, y) with A = (1 — S)/{1 + 6) and y = a/(l — 6). 
The latter parametrization makes explicit that for A = 0, the model consists of decoupled dimers: we take this to be 
our unperturbed Hamiltonian Hq. The rest of the Hamiltonian can be treated as a perturbation, and we can expand 
various physical quantities in powers of A. The formalism for studying n-particle sectors in perturbation theory is 
discussed in detail in a companion paper |L5| . 
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The series expansions for the ground state energy and triplet excitation spectrum have previously been computed 
[11 1 up to order 23. The two-particle excitations have been discussed using a leading order Brueckner ansatz calculation 
[13 1, a second order series expansion [O and an RPA study With our new technique, we perform high-order 

series expansions in powers of A for fixed values of y |ll|. As discussed in a companion paper |15|, we first calculate 
an effective Hamiltonian in the two-particle sector 

£ 2 (i,j;k,l)Hk,l|ff cff |i,j) , (3) 
and then calculate the irreducible two-particle matrix element 

A 2 (i, j; k, 1) = E 2 (i, j; k, 1) - E (S iM 5^ + <5i,i5j, k ) - Ai(i, k)<5jj - A x (i, l)<5 j)k 

-A 1 (j,k) ( 5 u -A 1 (j,l)^ k , (4) 

where 5 refers to a Kronecker delta function and Ai is the one-particle irreducible matrix element. 

Here we will only concentrate on the expansions for the following two lines in the parameter space: (i) a = 0, 
corresponding to nearest neighbor interaction only, and (ii) a = (1 — S)/2, which is a special line in the parameter 
space where the ground states are known exactly, also known as the Shastry-Sutherland line. The model at 5 = 
(A = 1) is the uniform Heisenberg chain in case (i) and the Majumdar- Ghosh model in case (ii). 



III. BOUND STATES WITH FORCED DIMERIZATION 



In this section we study the small-A regime, where our strong coupling expansions are convergent. Thus a simple 
truncation of the relevant power series expansions leads to highly accurate results. We discuss the number of bound 
states with different quantum numbers as well as their binding energies, the range of k-values where the bound states 
exist, the coherence length associated with the bound-pair as well as the singularities at the critical wavevector, where 
the binding energy goes to zero and the bound state merges into the continuum. We first consider the model with 
a = and then a = (1 — S)/2. 



1. Case a — 



For the case without the second neighbor interaction (a = 0), the series for the irreducible two-particle matrix 
element A 2 have been computed up to order 7 for singlet bound states, and to order 11 for triplet and quintet states 
p6[ . The reason why the singlet series is computed to only 7th order compared to 11th order for the triplet and 
quintet states is that the singlet has the same quantum numbers as the ground state. Thus a much more elaborate 
orthogonalization method is required to implement the cluster expansion jl5| . 

In this model, we find two singlet (Si and S2) and two triplet (Ti and T 2 ) bound states below the two-particle 
continuum, and two quintet antibound states (Q\ and Q 2 ) above the continuum. The existence of the second pair 
of bound states has not been reported by previous calculations, most likely due to a limited precision or a general 
incapability to deal with multiple bound states. The series for their energies (and also the lower edge and upper edge 
of the continuum) at band maximum k = n/2 are given in Table [i]. Note that there are some discrepancies for the 
energy of the lowest singlet bound state with the previous second order calculations |l3| . Our second order results 
agree with the series results of Barnes et al. jl4| , but disagree with the results of the Brueckner ansatz calculation [jl3) . 
Although the Brueckner ansatz is an expansion of the self-energy in terms of the density of excitations, it normally 
can recover the first few order of the series expansion in A exactly: it does give the correct second order result for the 
triplet bound state. Jl3| 

In the limit A — > 0, the formation of bound/antibound states Si, T\ and Q\ is well known, due simply to the 
interaction of two triplets on neighboring sites, and the wave function for Si, for example, is 

1^ (*)> ~E e 2 ^- +l /2) [4{j)tUu + 1) + tUMu + 1) 4(i)4o- + i)]|o > , (5) 
3 

where |0) is the ground state at A = consisting of nonoverlapping spin singlets on each dimer, and t' a (j) is a triplet 
creation operator which excites the singlet at j-th dimer into a triplet state with S z = a, (a — — 1, 0, 1). 

It is interesting also to look at the structure of these new bound states 5* 2 , T 2 and Q 2 in this limit. To compute 
their wavefunctions, one needs to diagonalize the second order effective Hamiltonian in the two-particle sector, which 
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can be reduced to an infinite dimensional symmetric tridiagonal matrix. Our calculations show that in this limit, 
the bound/antibound states S2, Ti and Q2 only exist at k = ir/2, and their wave functions are (we take S z = as 
example) 

(*o) = --= E e2tk3 E fntiv - n )*-iC? + ») 

j n=l 

+« t _ 1 (j - n)t\(j +n) - 4(j - n)4(j + n)]|0 > , (6) 
|Vr 2 (fc)> = -= E ^ E - n)*LxC3 +n)- t^U - + n)]\0 > , (7) 

j n=l 
1 °° 

|^Q 2 (fc)) = -= E ^ E - n)tU(j + n) 

j n=i 

+*LiCj - «)40' + ») + 24 U - n)tl(j +n)]|0 > , (8) 
where the amplitude /„ for two triplets sitting at j — n and j + n (i.e. separated by distance 2n) is 




/„ = < -V3(-2)-" for T 2 (9) 



Thus the formation of these new bound/antibound states is due to an effective interaction between triplets separated 
by an odd number of singlet dimers. It appears that Si, T\ and Q\ are fully "localized" states in this limit, with 
wavefunctions extending only across a single pair of dimers, whereas the states S2, T% and Qi are "extended", with 
the tail of the wavefunctions decreasing exponentially with distance. Our calculations show that at finite A and for the 
particular case k = 7r/2, the wavefunctions for S2, T% and Q2 still only involve triplets separated by an odd number 
of singlet dimers, while the wavefunctions for Si, 2\ and Q\ only involve triplets separated by an even number of 
singlet dimers. Thus the spectrum splits into two decoupled sectors at k = n/2. It would be interesting to explore 
the dynamics behind this phenomenon in more depth. 

With the wavefunction, one can also compute the coherence length L defined by 

L = pi| (10) 

where fd is the amplitude for two triplets separated by distance d. Note that the coherence length L defined here is 
measured in units of 2a, where a is the lattice spacing. At the limit A — * 0, L = 1 for S\ and T\ as expected, while 
for 52 and T 2 , L is 32/15 and 8/3, respectively. 

The two-particle excitation spectrum and the inverse of the coherence length 1/L versus momentum k for a rather 
large dimerization S = 0.6 are shown in Figs. [I] and ||. One can see that the singlet bound state Si exists for the whole 
range of momenta (its coherence length L is finite also for the whole range of momenta although the coherence length 
at k — is very large, about 6.7706), while other bound/antibound states exist only in a limited range of momenta 
k > k c . The "critical momentum" k c for a given bound state can be defined by the inverse of the coherence length 
1/L tending to zero or by vanishing binding energy. Technically, the first approach may give more reliable results on 
k c . The results for the critical momenta k c versus <5 are given in Fig. 0, where we can see that in the limit <5 — » 1, the 
S2, T2 and Q2 states exist only at k — n/2. We can also get the first few terms in the series expansion for k c for T\ 
and Qi states. Up to order A 3 , the dispersions for the T\ and Qi states are 

„ , T/1 rN „ 3 A 9 A 2 89 A 3 / A 3 A 2 105 A 3 \ 

Et, J +5) = 2 1 h h cos (2k) 

11 v ' 4 32 128 V 2 16 128 / v ' 

A 2 3 A 3 A 37A 3 cos(6/c) 4 



, cos(4 k) 1 — '- + (9(A 4 ) , (11a) 

4 16 J v ' 128 v ' ' v ' 

^ /t/, cn „ 3A A 2 21A 3 f\ 3A 2 37A 3 \ 
E ^ J{l + 6 ^ 2+ --T2 + U8 + {-2-^6 + U8) COS(2fc) 
A2 13A3 \ tA1 s 23A 3 cos(6A:) ^. x4 . 

T + ^r cos(4fc) + — i28^ + 0(A } - (llb) 
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FIG. 1. The excitation spectrum of the Ji — ,h — 5 chain with 8 = 0.6 and a — 0. Beside the two-particle continuum (gray 
shaded), there are two singlet bound states (Si and S2) and two triplet bound states (Ti and Tb) below the continuum, and 
two quintet antibound states (Qi and Q2) above the continuum. The inset enlarges the region near k = it/2 so we can see S2, 
T2 and Q2 below/above the continuum. 




FIG. 2. The inverse of the coherence length 1/L versus momentum k for two singlet (Si and S2), two triplet (Ti and Tb) 
and two quintet (Qi and Q2) bound/antibound states of the Ji — J 2 — 8 chain with 8 = 0.6 and a = 0. The inset enlarges the 
region near k = it/2. 
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FIG. 3. The critical S versus momentum k for singlet, triplet and quintet bound/antibound states of the Ji — J2 — S chain 
with a — 0. The dotted lines are the results of Eq. (113) . 



With this and the series for the 1-particle triplet excitation spectrum, one can get k c as 

f 2tt/3 + 5A/(4V3)- 757A7(192\/3) + 0(A 3 ) for T x 
c ~ \ 2tt/3 + %/3A/4 + 15V3A 2 /64 + 0(A 3 ) for Qi 

and in the limit k —> fc c , the behaviour of the binding energy is 

E b /JX = 4(fc - /c c ) 2 [3/16 - 11A/128 + 591A 2 /512 + 0(A 3 )] 

+8(k - k c ) 3 [V3/32 + 113a;/(256V3) - 25A 2 /(768%/3) + 0(A 3 )] + 0[(k - k c 



(12) 



(13) 



for T\, and 



E b /J\ = 4(k-k c ) 2 [ 



3 

16 



21 A 
128 
1 



39 A 2 



lV3(k - k c ) 3 [— + — - 



32 



256 
37 A 37 A 2 
256 ~ 1024 



0(A 3 )] 



0(A 3 )] + 0[(k - k c ) 4 



(14) 



for Qi. Here one can see that the "critical index" for E b in the limit k — > k c is 2, independent of the order of expansion, 
so one expects that this is exact. The results of Eq. ([l2]) are also shown in Fig. [| We can see that in the limit 5 — > 1, 
k c = 7r/3 for T\ and Qi, and as d decreases, k c for Qi increases, while k c for T\ firstly increases, then decreases. 

Actually S\ does not always exist in the whole range of momenta: it does not exist at k = when 5 — * 1. The 
inverse of the coherence length L for Si at k = versus <5 is given in Fig. [|, where we can see that as S approaches 
1, L diverges. Because A = is a critical point, we cannot get the series directly for the energy gap of S\ at k = 0. 
We can only get numerical results for it by solving the integral equation Jig ]. This makes it poorly convergent as 
A — > 1, as we will see in the next section. In this figure, we also plot the coherence length L at k — n/2 for Si and Tj 
(< = 1,2). < 

In the limit A — > (5 — > 1), the binding energy at k = n/2 for S\ and T\ is proportional to A, as expected, while 
for 5*2 and T2, the binding energy is proportional to A 2 . The rescaled binding energies E^/JX 1 versus S for Si and Tj 
(i = 1,2) are shown in Fig. ^| We also show some numerical exact diagonalization results ]l4| in this figure, which 
are in very good agreement with our series results. 

As evident from Figs. |i| and |^ the bound state T2 may disappear at about 8 = 0.2. 
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FIG. 4. The coherence length L versus 8 for two singlet (Si and 52) and two triplet (Ti and T2) bound states of the Ji — J2 — <5 
chain at fc = 0, 7r/2 and a = 0. The inset plots 1/L versus 8 for Si at k — 0. The results of the three highest orders are plotted. 
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FIG. 5. The scaled binding energy E^j JX n at k — tt/2 versus dimerization 8 for two singlet (Si and S2) and two triplet (Ti 
and T2) bound states of the Ji — J 2 — 8 chain with a = 0. The solid points are the numerical exact diagonalization results jL4j , 
The inset enlarges the region for T2. Several different integrated differential approximants to the series are shown. 
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FIG. 6. The energy gap E/J at k = tt/2 for two singlet (5*2 and S3) and three triplet (7i, T2 and T3) bound states versus 
the inverse of order 1/n for the Ji — J2 — S chain with 8 = 0.4 and a = (1 — 5)/2. The gray shaded regions are the two-particle 
continuum. The inset enlarges the region near k = tt/2 so we can see S3 and T3 below the continuum. 



2. Case a = (1 - 5)/2 

Along the special line a = (1 — <5)/2, the ground state is an exact product state, with the spins on each strongly- 
coupled bond forming a singlet. For non-zero S, the elementary excitations for this system are triplets. When 5 — > 0, 
the system has two degenerate ground states and the triplets unbind into a pair of free spin-half excitations. These 
spin-half objects, which are domain walls between the two ground states, are called solitons and they become the 
elementary excitations of the system. 

The series for the irreducible two-particle matrix element A2 has been computed up to order A 19 for 2-particle 
singlet, triplet and quintet states JlfjJljJ by using both orthogonal transformation (two block method) and similarity 
transformation methods fil| . Hence one can compute the series directly for the dispersion of the bound states using 
a degenerate perturbation expansion. It turns out that both transformations give identical series for the dispersion 
of the bound states, up to the order computed, although the series for the irreducible two-particle matrix elements 
A2 are different (the basis states are different in the two methods). The energy gap at k = tt/2 for one of the singlet 
bound states, Si, is 1 + 3(5 exactly The series for the energy gaps of the other bound states (and also the lower 

edge of the continuum) at k = and tt/2 are given in Table 0. 

Here we find three singlet (Si, S2 and S3) and three triplet, (T\, T2 and T3) bound states below the two-particle 
continuum, and two quintet antibound states (Qi and Q2) above the continuum. The dispersions for these bound 
states at S — 0.4 have been shown in Fig. 4 of a preceding Letter |fT9| . To demonstrate the reliability of our results, we 
plot in Fig. ^ the energy gap at S — 0.4 and k — n/2 for all bound states and the two-particle continuum versus the 
inverse of the order n up to n = 19. The results for Si are not plotted, since this case is known exactly. From this 
figure, we can see that the results are very well converged for n > 10. In the limit of A —> 0, the binding energy at 
k = tt/2 for Si and T\ is proportional to A, as expected, while for Si, Ti, S3 and T3, the binding energy is proportional 
to A 2 , A 2 , A 4 , A 6 respectively, as we can see from Table [n[ 

In the limit A —> 0, the wave functions for Si and Tj (i = 1, 2) are trivial, consisting of two triplets separated by 
i — 1 singlet dimers, while the wave functions at k — tt/2 for S3 and T3 are (here again we just take S z = for T3 as 
example) 

1 00 

\i, S M) = -iEE e2ife(i+n+1/2) /«[*tO')*-iO'+^+ 1 ) 
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FIG. 7. The amplitude /„ versus n for T 3 with 5 = 0.4, 0.6 and 0.9, and k = tt/2 



+tli(j)t\(j + 2n+ 1) - 40')40' + 2n + 1)]|0 > , (15) 

- oo 

|^(fc)) = -= £ £ ^ fc °' + " +1/2) /n[40> t - 1 (j + 2n + 1) - *Lx(i)4(j + 2« + 1)]|0 > . (16) 

For 5*3, the amplitude /„ for two triplets sitting at j and j ' + 2n + 1 is 

/„ = 2~ n V3 , (17) 

which decreases by a factor of 2 as n increased by 1. We cannot obtain an analytic expression for f n for T3, but the 
numerical results for 5 — 0.4, 0.6 and 0.9 are presented in Fig. [7| for small A (large S), f n is almost independent 
of n, so one has an infinite coherence length. From the above results, one can see that the bound states S3 and 
T3 are due to the effective attraction between triplets separated by an even number (exclude 0) of singlet dimers. 
Thus we find in this case that Si and Ti (i = 1,2) arc fully localized states, whereas S3 and T3 are "extended", with 
exponential tails to their wave functions. As for the case a = 0, our calculations again show that for k = tt/2 and 
any A, the wavefunctions for S2 and Ti only involve two triplets separated by an odd number of singlet dimers, while 
wavefunctions for Si, Ti, S3, and T3 only involve two triplets separated by an odd number of singlet dimers. 

The two-particle binding energies and the inverse of the coherence lengths 1/L for S = 0.4 are shown in Figs. || and 
^. One can see that the singlet bound state Si and the triplet bound state T\ exist for the whole range of momenta, 
while other bound states exist only in a limited range of momenta k > k c . The results for the critical momenta k 
versus 8 are given in Fig. |To|, where one can see that in the limit 6 — * 1, the S3 and T3 states exist only at k = tt/2. 
For S2 and Ti bound states, as before, one can get the first few terms in the series expansion for k c : 

_ / 2 1 /2a + 0.53033A 2 + O(A 3 ) for S 2 . . 



2tt/3 - 7A/(8V3) + 287A 2 /(768\/3) + <3(A 3 ) for T 2 

and in the limit k k c , the behaviour of the binding energy is 

E b /J\ 2 = 4(fc - A: C ) 2 [A 2 + 0(A 3 )]/32 

+8(k - fc c ) 3 ( — + 0.107723 A 2 + 0(A 3 ) ) + 0[{k - k c ) 4 } (19) 
V 32 y2 / 

for So, and 
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FIG. 8. The binding energy Et, for three singlet (Si, S2 and S3) and three triplet (Ti, Tb and T3) bound states versus 
momentum k for the Ji — J2 — 8 chain with 8 = 0.4 and a = (1 — <5)/2. The inset enlarges the region near k = 7r/2 so one can 
see the nonzero binding energy for S3, T3 . 
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FIG. 9. The inverse of the coherence length 1/L versus momentum k for three singlet (Si, S2 and S3) and three triplet (Ti, 
T2 and T3) bound states of the Ji — J2 — 8 chain with 8 — 0.4 and a = (1 — S)/2. The inset enlarges the region near k = 7r/2. 
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FIG. 10. The critical 8 versus momentum k for singlet and triplet bound states of the Ji — , J% — S chain with a = (1 — S)/2. 
The dotted lines are the results of Eq. (UM . 



E b /JX 2 = 4(fc - k c f 



+8(k-k c ) 3 



0(A 3 ) 



3 97 A 461 A 2 
32 + 512 + 1024 

/ V3 199 A 19273 A 2 
I 64 1024 %/3 24576 V3 



+ 0(A 3 ) +0[(k-k c ) 4 



(20) 



for T2. Here again one can see that the "critical index" is 2. The results of Eq. fllg| ) are also shown in Fig. 10. We can 
see that in the limit S — ► 1, k c = and 7r/3 for S2 and T2 respectively, and as <5 decreases, k c for S2 firstly increases, 
then decreases back to at around 8 — 0.38, while k c for T2 firstly decreases, then increases. We also can see from 
this figure that S3 and T3 only exist over a tiny range of momenta for all 6. 

As for Si in the case a — 0, 5 — * 1 is also a critical point for T3 at k = ir/2. The inverse of the coherence length 
L for T3 at k — n/2 and S2 at k — versus S are given in Fig. O, where we can see that as 5 approaches 1, the 
coherence length for T3 at k = ir/2 diverges, while the bound state S2 at k = appears at S < 0.38, consistent with 
Fig. Eol In this figure, we also plot the coherence length for other bound states at k = and ir/2. L is exactly 1 for 
Si at fc = 7r/2, while for T±, S2 and T2 at A: = tt/2, L is almost equal to 1,2,2 respectively for all S. For Si and T\ at 
k = 0, L is 1 in the limit 5 — > 1, and as 5 decreases, L increases. For S3 at fc = 7r/2, L is 11/3 in the limit S — > 1, and 
as <5 decreases, L decreases. 



IV. REGIME OF VANISHING BOND-ALTERNATION: UNBINDING OF SPIN-HALF EXCITATIONS 



In this section we turn to the regime of small S or A near unity. In this case, our results are less accurate and we 
have to rely on series extrapolation methods. The limit <5 — > is a critical point, where we expect singularities in 
various physical quantities. Hence, the convergence of the series breaks down and a simple truncation does not lead 
to meaningful results. In this case, we use the Dlog Pade and integrated differential approximants Q to extrapolate 
the series for the single-particle energies and the two-particle binding energies. We present results based on these 
extrapolations. 

We begin this section by making a few comments about the extrapolation of one and two-particle energies to the 
5^0 limit for the nearest neighbor model (a = 0). In this case the uniform system at 5 = is the Bethe-ansatz 
solvable nearest-neighbor Heisenberg model, with no gap in the excitation spectrum. Furthermore, it is believed 
that the mapping on to the massive Thirring model gives the exact spectrum at small k for small S pJ]. The latter 
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FIG. 11. The coherence length L at k — and tt/2 versus S for three singlet (Si, 5*2 and 5*3) and three triplet (Ti, T2 and 
T3) bound states of the Ji — J2 — S chain with a = (1 — 5)/2. The inset plots l/L versus 8 for S2 at fc = and T3 at fc = 7r/2. 
The results of 17th to 19th order series are plotted. 



model has a well defined singlet and a well defined triplet excitation whose energies are in the ratio of \fi. There 
are no further bound states so that there are no triplet states between the single particle gap A and the two-particle 
continuum gap 2A. 

As evident from figures ||, and || only one singlet bound state exists as 8 goes to zero for fc = 0. Its coherence 
length appears large but finite in our calculations, though its binding energy goes to zero. This, together with the 
elementary triplet whose energy also goes to zero as 8 goes to zero, gives the two well defined states expected from 
the massive Thirring model. As discussed in the previous section, the convergence of the singlet excitation energy at 
k = becomes poor as A = 1, as it has to be gotten from a numerical solution of the integral equation rather than 
from an extrapolation of a binding energy series. This ratio is plotted in Fig. |l2|. It clearly stays close to 2 down 
to (5 = 0.2 and begins decreasing for smaller values of 8. However, for 8 < 0.1 the extrapolations become completely 
unreliable and we cannot tell if the ratio approaches as 8 — > 0. 

We now turn to the frustrated model with a = (1 — S)/2. First, we show the dispersions of the single-particle 
triplet excitation, the lowest-energy two-particle singlet bound state Si, and the bottom of the two-triplet continuum 
over the full Brillouin zone for various values of 8 in Fig. [l3[ It is evident that the triplet and the singlet spectra 
become degenerate as 8 goes to zero. This is direct evidence for free spin-half excitations, since a pair of free spin-half 
excitations will form singlet and triplet states of equal energy. Note that previous series expansion studies have shown 
JTT| that the spectral weight associated with the triplets does not vanish over the entire Brillouin zone. It remains 
finite in the vicinity of fc = tt/2 at <5 = 0. This result is consistent with the variational calculations of Shastry and 
Sutherland, H who had found a soliton-antisoliton bound state near fc = tt/2. Shastry and Sutherland had also found 
that the bound states were four-fold degenerate meaning that there are degenerate singlet and triplet bound states. 
Thus, our results are completely consistent with their calculations. S0rensen and collaborators Q have also given an 
exact demonstration of triplet and singlet states with energies (1 + 8) J and (1 + 38) J at k = tt/2 respectively, and 
thus a splitting of 28 J which vanishes as 8 — > 0. Our results fit this pattern precisely. 

Note further from this figure that as 8 goes to zero, in the region not too far from fc = tt/2, the two-triplet continuum 
falls below the states Si and Ei (where Ei is the single-particle triplet excitation). This is also true for the variational 
calculation of Shastry and Sutherland shown in Fig. |ll](d), though, to our knowledge, it has not been noted before. 
At fc = tt/2, stable singlet and triplet states are rigorously known to exist. A plausible picture is that a stable state 
exists only at or very near this wavevector, and even at this point its binding energy with respect to the multiparticle 
continuum is extremely small. It is also likely that the spectral weight of these states is only appreciable in this narrow 
region; further from fc = tt/2, these states will be lost in the continuum, and their spectral weight will be neglible, in 
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FIG. 12. Ratio of the lowest singlet excitation energy Esi to the lowest triplet excitation Ei, in the nearest neighbor model 
as a function of 8. Results at order 4 (dotted line), 5 (short dashed line), 6 (long dashed line) and 7 (solid line) are plotted. 



agreement with the calculation of Singh and Zheng 

In the series extrapolation for the figure above, we have made use of the fact that as 8 goes to zero (or A goes to 
1), the energy gap approaches a constant with correction proportional to (1 — A) 2 / 3 ||, so in the series extrapolation 
we transform the series to a new variable 



A' = 1 - (1 - A) 



2/3 



(21) 



to remove the singularity at A = 1 . Byrnes et al. Q also predict that the singlet-triplet splitting in the limit 5 — > is 

{E Sl -E 1 )/J = c 1 S + c 2 8 5 / 3 . (22) 

Our series analysis seem to favor these arguments, as can be seen from Fig. |lj. But our results give c\ — 2.25(5), 
rather than 15/8 as given by Byrnes et al. this difference is probably because we choose a = (1 — S)/2, while in 
the calculations by Byrnes et al. a is fixed to be 1/2. 

In Fig. |l5] the binding energy for the lowest triplet bound state Ti is shown. It is clear that the extrapolation for T\ 
does not converge very well. Looking at the figure, it is plausible to suggest that as S goes to zero the triplet binding 
energy also approaches A ss 0.24J, the single particle energy gap. If true, this implies that this state also becomes 
degenerate with the single particle state in this limit. This is consistent with the idea that there are an infinite number 
of states between the single particle gap A and the two-particle continuum, which begins at 2A. Thus the lowest of 
these states must approach A in energy. However, we do not see a large number of two-particle bound states in our 
calculations. We suspect that the other bound states may arise from 3, 4, 5, . . ., particle states. We hope to extend 
our methods to study such multiparticle bound states in the future. 

We now turn to the bound states at k = ir/2. We show the rescaled singlet and triplet binding energies at k = ir/2 
in Figs. [l(] and [l?]. These results for Si (i = 1,2,3) and T,; (i = 1,2) are obtained from the integrated differential 
approximants to the series, while the results for T3 are obtained from the numerical solution of the integral equation 
J15I at orders 14 to 19, since we cannot get the series directly for T3. We can see from these figures that as S — ► 0, 
the binding energies for all Si (i — 1,2,3) approach the same value, which is close to or equal to zero. Among Tj 
(i = 1,2,3), T\ has the largest binding energy in the limit 6 —> 1, but at 6 = 0.092(3), the binding energy for Ti 
becomes zero, while the binding energies for Ti and T3 are still nonzero, i.e., there appears to be a level crossing 
between Tj (i = 1, 2, 3) here: the level crossing between T\ and T2 happens at S = 0.221(1) where the binding energy 
is Eb/ J = 0.11955(10). The reason that Ti and T2 can cross smoothly is presumably that the bound state for T\ only 
involves two triplets separated by an even number of singlet dimers, while the bound state for T2 only involves two 
triplets separated by an odd number of singlet dimers, so that the two states lie in disjoint sectors. 
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FIG. 13. Dispersion of the single-particle triplet excitation Ei (solid line), the lowest two-particle singlet bound state Si 
(dashed line), and the bottom of the two-triplet continuua (dotted line) for 8 — 0.1 (a), 0.02 (b) and (c) for the Shas- 
try-Sutherland model. Window (d) gives the variational results of Shastry and Sutherland, where only the continuua are 
shown. 
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FIG. 14. The 1-particle triplet gap (E\) and lowest singlet energy gap Es 1 at k — (upper window), and (Es 1 — E\)/J5 
(lower window) versus 8 2 ^ 3 for the Ji — J 2 — S chain with a = (1 — <S)/2. The results of several different integrated differential 
approximants to the series are shown. 



0.3 




FIG. 15. The binding energy E b /J of the lowest triplet bound state (Ti) at k = as a function of 5. The results of several 
different integrated differential approximants to the series are shown. 
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FIG. 16. The rescaled binding energy Eb/J\ n at k = tv/2 versus dimerization <5 for three singlet bound states Si (i = 1, 2, 3) 
of the Ji — J 2 — S chain with a — (1 — 8)/ 2. The inset enlarges the region near 5 = 0. Several different integrated differential 
approximants to the series are shown. 




FIG. 17. The rescaled binding energy Eb/JX n at k = n/2 versus dimerization S for three triplet bound states Ti (i = 1, 2, 3) 
of the Ji — J 2 — 5 chain with a = (1 — S)/2. The inset enlarges the region for T3. For Ti and T2, the results of several different 
integrated differential approximants to the series are shown, while for T3 the numerical results from the integral equation at 
order 14 to 19 are shown. 
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V. CONCLUSION AND DISCUSSION 



In this paper we have carried out an extensive investigation of the two-particle spectra of the frustrated alternating 
Heisenberg chains using strong coupling series expansions. 

In the regime of weakly coupled dimers, the elementary excitations of the system are triplets and the spin-half 
excitations are confined. In this regime our series expansions are convergent, and we have studied in great detail the 
properties of the two-particle spectra, including binding energies, coherence lengths and critical properties associated 
with vanishing binding energies and diverging coherence lengths. We find in every case, just as for the spin ladder 
system ]l5| , that where a 2-particle bound state emerges below the continuum at a "critical momentum" k c , the 
binding energy behaves like (fc — k c ) 2 as k — > fc c+ , and the coherence length diverges, as one would expect. Several 
distinct bound states can be identified, particularly near k = tt/2. Many of the bound states can only be seen by 
going to sufficiently high orders in the perturbation expansion, showing the extended character of the pair-attraction. 

We have also studied the regime of deconfined spin-half excitations (5 — ► 0) by using series extrapolation methods. 
Several properties of one and two-particle spectra give a clear indication of this deconfinement transition. The spectral 
weights of the triplets vanish at the transition (except near k = 7r/2) and the singlet and triplet excitations become 
degenerate. These methods can be used to look for such unbinding transitions in higher dimensional models as well. 

Our studies also raise several puzzles that need to be addressed in the future. How does the spectrum between A 
and 2A get filled up as one goes from the confined to the deconfined phase? We suspect that multi-particle bound 
states, with varying number of triplets are important. This needs to be further addressed. The picture is much clearer 
in the soliton language, as one might expect: at small S, there is a discrete spectrum of s — s bound states confined 
by a linear potential, which becomes continuous as S — > and the confining potential vanishes Another puzzle 

is the crossing of energy levels. As A goes from zero to unity along the Shastry-Sutherland line, one expects to see 
crossing of n-particle states, with different n. How does this take place? How can this be accounted for within the 
perturbation theory? Finally, are the Shastry-Sutherland bound states stable away from k — tt/2, where they might 
decay into the 4-soliton continuum? We hope to address some of these issues in future. 
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TABLE I. Series coefficients for dimer expansions of the energy gap E/J(l + 5) of two singlet bound states (Si and S2), 
two triplet bound states (Ti and T 2 ), two quintet antibound states (Qi and Q2), and the lower edge and upper edge of the 
continuum (C; and C u ) at k = tv/2 for the Ji — J2 — 5 chain with a = 0. Nonzero coefficients A" up to order n = 11 are listed. 



n 


E S , /J(l + S) for Si 


Es 9 JJ(l + S) for S 2 


E T JJ(1 + S) for Ti 


E T9 /J(1 + 8) for T 2 





2.000000000 


2.000000000 


2.000000000 


2.000000000 


1 


-5.000000000X10" 1 


0.000000000 


-2.500000000 xlO -1 


0.000000000 


2 


1. 875000000x10^ 


-3. 906250000x10^ 


1.562500000 xl0 _1 


-2.812500000X10 -1 


3 


3.906250000xl0~ 2 


7.812500000xl0~ 2 


-2.343750000xl0~ 2 


6.250000000xl0~ 2 


4 


-2.050781250xl0~ 2 


3.324584961 xl0 _1 


-6.168619792xl0~ 2 


9.326171875xl0~ 2 


5 


-4. 154459635 xl0~ 2 


-5.464426676 xl0~ 2 


-8.189561632xl0~ 2 


-1.143256293xl0~ 1 


6 


-4. 004075792 xl0~ 2 


-4.499045478 X 10" 1 


-9.442308214xl0~ 2 


-1.716450585X10 -1 


7 


-2.573965214xl0~ 2 


-1. 200470705x10^ 


-8.740929615 xl0~ 2 


-2.332303701 xl0~ 2 


8 






-6.322668408xl0~ 2 


2.681890517xl0~ 2 


9 






-2.250317901xl0~ 2 


-3.523763254xl0~ 2 


10 






2.941323970xl0~ 2 


1.450322827xl0~ 2 


11 






9.059538250X10 -2 


1.467690427xl0~ 1 


n 


E Q2 /J(1 + S) for Q 2 


E Ql /J{l + 5) for Qi 


E C ,/J(1 + S) for Ci 


E Cu /J(l + 5) for C u 





2.000000000 


2.000000000 


2.000000000 


2.000000000 


1 


0.000000000 


2.500000000 xlO -1 


0.000000000 


0.000000000 


2 


3.125000000X10 -2 


-9.375000000xl0~ 2 


-2.500000000 XlO -1 


0.000000000 


3 


1.250000000X10- 1 


-1.015625000X10- 1 


3.125000000xl0~ 2 


1.562500000X10" 1 


4 


1.513671875X10- 2 


-5.126953125xl0~ 2 


-2.343750000xl0~ 2 


9.375000000xl0~ 2 


5 


1.916503906X10- 2 


-1.114908854X10- 2 


-6.868489583xl0~ 2 


-7.269965278xl0~ 3 


6 


5.786874559xl0~ 2 


-6.629096137xl0~ 4 


4.757351345 xl0~ 3 


-1.196628147xl0~ 2 


7 


2.270891195xl0~ 2 


1.044698998xl0~ 2 


-1. 730789373 xl0~ 3 


8.809831407X10" 3 


8 


1.947733239xl0~ 2 


2.007299524xl0~ 2 


-3.478056560xl0~ 2 


4.201953303 xl0~ 3 


9 


2.277670953xl0~ 2 


3.177196538X10" 2 


2.936064788xl0~ 4 


-2.457848151 xl0~ 3 


10 


-1.225328947xl0~ 3 


4.251490921 xl0~ 2 


5.560815615xl0~ 3 


8.842688260X10- 4 


11 


-9.435317910X10- 3 


5.454945223 xl0~ 2 


-2.004222569xl0~ 2 


2.355556912xl0~ 3 
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TABLE II. 


Series coefficients for dimer expansions of the energy g 


ap E/J{\ + 5) of singlet bound states (Si, S2, S3), triplet 


bound states (Ti and T2) and the lower 


edge of the continuum (C;) at 


k = and ir/2 for the Ji - J 2 


-5 chain with a = (l-<5)/2. 


Nonzero coefficients of X" up to order 


n = 19 are listed. 






n 


Si at k = 


Ti at k - 


Ci at k - 


S 2 at k = 7I-/2 





2.000000000 


2.000000000 


2.000000000 


2.000000000 


1 


-1.000000000 


-5.000000000X10" 1 


0.000000000 


0.000000000 


2 


-5.000000000X10 -1 


1.250000000X10" 1 


-1.000000000 


-7. 500000000x10^ 


3 


-2. 500000000x10^ 


-1. 562500000x10^ 


-5.000000000X10" 1 


-6.250000000X10" 1 


4 


6.250000000xl0~ 2 


-9.257812500X10- 1 


6.250000000X10- 2 


-4.687500000X10- 2 


5 


2.031250000xl0~ 1 


-1.833007813 


3.593750000X10- 1 


8.066406250X10- 1 


6 


3.255208333xl0~ 2 


-2.977783203 


1.542968750X10- 1 


1.142822266 


7 


-2.444118924xl0~ 1 


-4.028717041 


-2.992621528X10 -1 


-2.324761285X10 -2 


8 


-2.273898655X10 -1 


-3.409357212 


-4.475063748X10" 1 


-2.524002923 


9 


1.677377960xl0~ 1 


4.411956160 


-2.119700114xl0~ 2 


-3.662303224 


10 


4.367628576X10" 1 


3.284248577X10 1 


5.406151312X10" 1 


6.600969298X10" 1 


11 


4.447491734X 10~ 2 


1.064643399X10 2 


4.678734733X10- 1 


1.029881729X10 1 


12 


-6.649650688X10" 1 


2. 537572572 x 10 2 


-3.650978943X10" 1 


1.420401469X10 1 


13 


-5.647530897xl0~ 1 


4.601292420X10 2 


-9.976553176X10 -1 


-5.430672245 


14 


6.692916069xl0~ 1 


5.089128172X10 2 


-3.568290713X10 -1 


-4. 753274080 xlO 1 


15 


1.434455637 


-4.198436299X10 2 


1.149565251 


-6.050262051 xlO 1 


16 


-7.932368314xl0~ 2 


-4. 530966202 x 10 3 


1.548592323 


3.784212529X10 1 


17 


-2.535956672 


-1.638065807X10 4 


-3.749497848xl0~ 1 


2.350096418X10 2 


18 


-1.777648180 


-4.239290399X10 4 


-2.734332838 


2.710161551X10 2 


19 


3.125549537 


-8.305280519X10 4 


-1.897635556 


-2.494277528X10 2 


n 


S3 at k = 7r/2 


Ti at k - ir/2 


T 2 at k - k/2 


Ci at k = tt/2 





2.000000000 


2.000000000 


2.000000000 


2.000000000 


1 


0.000000000 


-5.000000000X10- 1 


0.000000000 


0.000000000 


2 


-5.000000000X10- 1 


-1.250000000X10- 1 


-6.250000000X10- 1 


-5.000000000X10- 1 


3 


-2. 500000000x10^ 


0.000000000 


-4.375000000X10" 1 


-2. 500000000x10^ 


4 


-7.421875000xl0~ 2 


-3.906250000xl0~ 3 


-2.148437500X10- 1 


-6.250000000xl0~ 2 


5 


1.123046875xl0~ 2 


-6.835937500xl0~ 3 


2. 929687500 xl0~ 3 


3.906250000xl0~ 2 


6 


5.106608073xl0~ 3 


-5.777994792xl0~ 3 


1.789143880X10" 1 


4.231770833xl0~ 2 


7 


-4.116566976X10- 2 


-8.443196615xl0~ 4 


3.243408203X10- 1 


-3.363715278xl0~ 3 


8 


-6.972208729xl0~ 2 


-4.641285649X10" 3 


4.088701319X10- 1 


-4.214590567xl0~ 2 


9 


-5.548241367xl0~ 2 


2.439375277xl0~ 4 


3.233009297X10" 1 


-4.348302771 xlO" 2 


10 


-1.615796414xl0~ 2 


-5.090890892xl0~ 3 


-7.412812403xl0~ 2 


-1.641364313xl0~ 2 


11 


1.393358333xl0~ 2 


2.287079761 xl0~ 3 


-7.776608991 xlO" 1 


9.839517861 xlO" 3 


12 


1.628315511X10" 2 


-6.726104563xl0~ 3 


-1.459665977 


1.402656002xl0~ 2 


13 


1.092421645xl0~ 3 


4.895928063xl0~ 3 


-1.501800181 


-1.566731983xl0~ 3 


14 


-7.495980128 xl0~ 3 


-9.649045824xl0~ 3 


-3.338269791 xlO" 1 


-1.834857746xl0~ 2 


15 


4. 106539880 xl0~ 3 


8.959512690xl0~ 3 


2.080412782 


-2.058150001 xlO -2 


16 


2.756737781 xl0~ 2 


-1.466894655xl0~ 2 


4.898068211 


-8.451718109xl0~ 3 


17 


4.306603014xl0~ 2 


1.550047623xl0~ 2 


6.361874618 


4. 993840865 xlO -3 


18 


3.910407449xl0~ 2 


-2.303464950xl0~ 2 


4.129539836 


7.616721351xl0~ 3 


19 


2.231265079xl0~ 2 


2.621842144xl0~ 2 


-3.760012157 


-1.055977079xl0~ 3 
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